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Abstract 



A relativistic covariant one boson exchange model, previously applied to 
describe elastic nucleon-nucleon scattering, is extended to study r] production 
in NN collisions. The transition amplitudes for the elementary BN — ► 7]N 
processes with B being the meson exchanged ( B = ir,(T,r], p,u> and S) are taken 
to be the sum of four terms corresponding to s and u-channels with a nucleon or 
a nucleon isobar N* (1535 MeV) in the intermediate states. Taking the relative 
phases of the various exchange amplitudes to be +1, the model reproduces the 
cross sections for the NN — > Xt] reactions in a consistent manner. In the limit 
where all r/'s are produced via N* excitations, interference terms between the 
overall contribution from the exchange of pseudoscalar and scalar mesons with 
that of vector mesons cancel out. Consequently, much of the ambiguities in 
the model predictions due to unknown relative phases of different vector and 
pseudoscalar exchanges are strongly reduced. 
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1 Introduction 



Much interest, both experimental [1-10] and theoretical [11-26], was devoted in recent 
years to the production of rj meson in hadronic collisions and electroproduction pro- 
cesses. In analogy with rr production through excitations of the A (1232 MeV) P33 
resonance, it is believed that 77-meson production in nucleon-nucleon (NN) collisions 
proceeds via the excitation of nucleon isobars. Because of isospin conservation how- 
ever, only I=| N* isobars which decay into 77N pair are important. Pion and photon 
induced reactions^], [27| provide strong evidence that near threshold, the N*(1535 
MeV) Sn resonance dominates 77 production. This is not unexpected in view of the 
proximity of its mass to that of the 77N pair and, its large branching ratio ( 30-55% 
) for N* — > rjN decay p8|. The couplings to other resonances are week and seem to 
play a minor role close to threshold. 

The production of 77 mesons in NN collisions was considered by several groups [11- 
16] within one boson exchange (OBE) models, where the N* (1535 MeV) is excited via 
the exchange of a boson B and then decays into 77N pair. Much of the uncertainties 
in the predictions of these model calculations concern the coupling of the mesons 
exchanged to the Sll resonance. The coupling of the 7r and 77 are known from the 
decay widths of this resonance but not so well the coupling of other mesons. Assuming 
vector meson dominance, Germond and Wilkin|12[ determined the p coupling from 
7iV — > N* data and, as there is no evidence for isoscalar photon coupling to the N* 
they neglect the contribution from u exchange. Vetter et al. ||14|| determined the p 
coupling from the empirical decay width of 5% for the N* to decay into two pions. 
They suggested a value q p nn* = 0.615 , nearly a factor of ~ 2.7 weaker than that 



deduced by Germond and Wilkin |T2[. Although no experimental evidence exists, they 
assumed that the uj also couples to the N* with g^NN* = 0.436, resulting with about 
equal p and u contributions]!!]]. Contrary to this, with the p coupling taken from 
vector dominance, the cross section is dominated by p exchange ||16|| while the u with 



a coupling as proposed by Vetter et al.|Tj] may influence the total cross section at 
most by ~ 30 — 40%. It is not clear though whether such a contribution is required 
to explain the existing pp — > ppr] data. 
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In addition to these uncertainties concerning the coupling constants, the relative 
phases of different meson exchange amplitudes are unknown. Taking these to be ones 
of signs and restricting the model to tt, rj, uj and p exchanges only, there are eight 
solutions for the transition amplitude, one for each of the sign combinations possible. 
A solution, where the p and n exchanges add destructively, seems to explain best the 
cross section data for the pp —>■ pprj reaction jjTjl O , [IB]. 

It is the aim of the present work to consider rj meson production in NN collisions 
within a relativistic covariant OBE model, where a virtual boson B produced on one 
of the incoming nucleons, is converted into an 77 meson on the second nucleon (see 
Fig. 1). Each of the different exchange contributions to the production amplitude is 
determined by three factors, representing a source function of the B-N-N vertex, a 
boson propagator and a conversion amplitude for the BN — ► r]N process. Only the 
latter depends on unknown boson-isobar couplings. 

In variation with previous calculations the 77 production amplitude includes four 
terms, corresponding to s and u- channels with nucleon isobars (diagram 2a-2b) and 
nucleon excitations (diagram 2c-2d) in the intermediate states. At threshold the con- 
tribution from s-channel N* pole (diagram 2c), hereinafter referred to as a resonance 
production term, dominates the process. The other diagrams, by far less important, 
furnish a background term. Though the various exchange amplitudes and calculated 
cross sections remain similar to the ones calculated previously using a non relativistic 



OBE model|16[, the present model has distinct features which make predictions con- 
siderably more reliable. Particularly, in the limit where all 77's are produced via N* 
excitations, interference terms between the overall contribution from scalar and pseu- 
doscalar meson exchanges with that from vector mesons cancel out. Consequently, 
the model predictions are insensitive to the choice of the relative phases of vector 
meson exchanges with respect to those of scalar and pseudoscalar mesons. 

The paper is organized as follows. In Sect. 2 details of the formalism and model 
parameters are presented. In Sect. 3 we write the S-wave production amplitude in 
a simple form suitable for numerical calculations. In Sect. 4 we discuss interference 
between various exchange contributions. The relative importance of various contri- 
butions to the amplitude and cross section for the NN — > NNrj reactions is studied 
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in Sect. 5. At energies close to threshold, final state interactions (FSI) strongly influ- 
ence the scale and energy dependence of the calculated cross sections[[L5], [16| . These 



are introduced in Sect. 6 where comparison with data is to be made. We conclude 
and summarize in Sect. 7. 

2 Formalism 

To calculate the primary production amplitude for the NN — > Xt] reactions, with 
X designating a bound or unbound two-nucleon state, we employ a fully covariant 
formalism based on an effective OBE model. We assume a reaction mechanism as 
depicted in Fig. 2, where a boson B created on one of the incoming nucleons is 
converted into an r\ meson on the second one via a conversion process, BN — > r)N. 
We include both nucleon and N* (1535 MeV) isobar excitations in the intermediate 
states. We do not include diagrams in which the 7] is produced on an internal meson 
line. Because of isospin and parity conservation only diagrams with a and S lines may 
contribute, but these are expected to be negligible small due to the weak coupling 
constants involved. 

2.1 Lagrangian 

The effective Lagrangian density is taken as, 

L = l^N^dJrN + ln^LN^duVN + g a NN^a N + 

g S NNNr5N + N[g pNNl » + l^a^d v ]rp,N + g uNN Nj^N + 
[g^NN* N*T7rN + g V NN* N*r]N + ig aNN ,N*^ 5 aN + 

g P NN*N*l 5 Y?P,N + g^ NN *N*^Yu,N + ig SNN *N*^5N + h.c] , (1) 

where N, N*, n, a, rj, p, uj and 5 represent the fields of a nucleon, a nucleon isobar and 
the mesons in the spin isospin space; M and itlb are the the mass of the nucleon 
and a meson B; r are Pauli matrices acting in the isospin space and the 7's denote 
Dirac matrices. This is a simple generalization of a Lagrangian used by Machleidt et 
al.p^j to describe NN elastic scattering data. As in Ref.p9| a pseudovector coupling 



is assumed for the 7rNN and 7/NN vertices and the tuNN tensor coupling f^NN is set to 



be zero. It is to be noted that using pseudoscalar couplings rather than pseudovector 
couplings for the 7r and rj mesons would not alter the principal conclusions from the 
present work. But, taking the u tensor coupling to be zero has the consequence that 
the uj exchange contribution to rj production rate via N* excitations vanishes for the 
pp — > pprj reaction (see Sect. 3). 

Note also that the interaction Lagrangian densities for the a and 5 mesons in Eqn. 
HI are chosen in the form igsNN* [N*~f 5 NB+N~f 5 N*B]. As we show in Appendix A, the 
coupling constants (Jbnn* (B — a,S) are related through an effective triangle diagram, 
to QbnN) g-wNN and g n NN*, and this form ascertains that all of these quantities be real 
numbers. 

2.2 The primary production amplitudes 

Suppose now that the production of a pseudoscalar meson NN — > NNP reaction 
proceeds via the mechanism depicted in Fig. 1. The primary production amplitude 
can then be written in the form, 

M$ = T,{ T BN^PN(p4,k ] p 2 ,q)G B (q)S BNN (p 3 ,p l )] + [1 <-> 2;3 <-> 4], (2) 

B 

where Sbnn{P3,Pi) and Gb(q) represent the source function and propagator of the 
boson B, and Tbn-*pn being the amplitude for the BN —>■ PN transition. For 
the reactions to be considered, the latter depends on the unknown couplings to the 
nucleon isobar. The parameterization of the source functions and boson propagators 
is rather well determined from fitting NN elastic scattering data. In Eqn. ^Pi,q and 
k designate 4-momenta for the i-th nucleon, the boson exchanged (B) and the meson 
produced (P). The sum runs over all possible boson exchanges which contribute to 
the process. The bracket [1 *-> 2; 3 «-> 4] stands for a similar sum with the momenta 
Pi,Ps and P2,P^ being interchanged. 

To make the present paper reasonably self contained we write in what follows 
the expressions for the source functions, propagators and conversion amplitudes. We 
use covariant parameterizations for Sbnn and T BN ^ PN . For scalar and pseudoscalar 



mesons the propagator and source function are[pQ, 31 



G s Xq 2 )=i/(q 2 -rn 2 B + ie) , (3) 



Ssnn = u(p 3 )Iu( Pl )F s (q 2 ) , (4) 
Spnn = u{pz)^Iu{pi)F P (q 2 ) , (5) 

where Fp{q 2 ) is a source form factor and I the appropriate isospin operator, i.e., 
1=1 and I = t for isoscalar and isovector mesons, respectively. Here u is a nucleon 
Dirac spinor and p^ = p\ — q is the final nucleon momentum (see Fig. 1). The Dirac 
spinors as well as the amplitudes T BN ^ NN and M?j? are normalized as in Itzykson 



and Zuber|3l|| 



Without making unnecessary limiting assumptions about the reaction mechanism, 
the amplitude for the conversion process of a scalar into a pseudoscalar meson , 



SN — > PN, has the usual form [BO 



T sn ^pn(P4, k;p 2 , q) = u{p A )^[A s H —B s ]Iu(p 2 ) ■ (6) 

Similarly, the amplitude for the conversion of a pseudoscalar meson into another 
pseudoscalar meson, P X N — > P 2 N, is given by, 

T PlN ->p 2N (p4,k;p2,q) = u(p 4 )[A P + -(ft + $)B P ]Iu(p 2 ) . (7) 

Here ^ = and the quantities A$(p), B S (p) are invariant functions depending on 
the Mandelstam variables. These functions have an isospin structure in accordance 
with the following rules : 

i) As = A x 1 for both B and P being isoscalar particles; 

ii) As = A x t for B being isovector and P isoscalar or vice versa; 

iu)A s b = A x 5 ab + A 1 x [r a ,r b ] for B and P being both isovector particles. In 
Subsection. 2.3 we shall calculate these functions for the model proposed using the 
effective interaction Lagrangian, Eqn. |1| 

The propagator for a vector meson is defined as, 



G 



q 2 — m 2 B + ie 



where g^ v is the metric tensor. Generally, a vector source is the sum of vector and 
tensor current terms, 



S% n (Pi,Ps) = u{p,)[rF$\q 2 ) + ia^q u F$\q 2 )]Iu{p 2 ) , (9) 
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with Fy\q 2 ) and Fy\q 2 ) standing for vector and tensor source form factors, quan- 
tities being the analogous of the nucleon electromagnetic form factors. This last 
expression contains conserved currents only and therefore satisfies current conserva- 
tion, i.e., 

q^S v NN{q) = . (10) 

Likewise, the conversion amplitude of a vector meson into a pseudoscalar meson, 
VN — > PN, is expressed in terms of eight invariant functions Ay ti (i = 1 — 8) all 
having an isospin structure as mentioned above for the As and Bs- More concisely 
we write, 

T%^ PN {p 4 ,k;p 2 ,q) = u( Pi ) 1 5 {YAV 2+ p%A^ + k»AV 6 + q»AV 8 ]Iu( P2 ) , (11) 

where A\- = Ay t i — ftAyj. It is worth noting that because of current conservation, 
Eqn. [TU|, only the first term of Eqn. 8 may contribute to the production process so 
that, 

^liu^VNN - — 2 ™2 T Z VNN ' ^ > 

q "'k t tt 

Also, for this same reason, the term q^A^ 8 in Eqn. |TT] does not contribute to the 
production amplitude, M^ n \ By using straightforward algebra of 7 matrices and the 
Dirac equation for a free particle and, by substituting equivalent two-component free 
spinor matrix elements (Table A6.1 of Ref. |3^]) in Eqns. [3|-|TTl we ma Y write the 
primary production amplitude as, 

Mp N ^ PNN (pi,P2\ Pa, P4, fc) = £ M s n + Yl M p + X^ n + [1^2;3^4], (13) 

S P V 

where the sums over S, P, and V run, respectively, over all scalar, pseudoscalar and 
vector mesons and, 

M™ = S + S k 13 a k 3 + S* + S kl a k ,o l 2 , , (14) 

MP = P&U + p u 4A , (15) 
My =V + V k 3 a k 13 + V 2 \a k 2 , + V H a k 3 a l M . (16) 



In Eqns. p^^ - flq all of the quantities S ,V ; S^-, P^, V^-; Sy, V^; S r ' 
functions of the particle momenta involved. Their definitions are given in Appendix 
B for the general case in terms of the invariant functions. 



2.3 The invariant functions 

To complete the derivation of the i] production amplitude the invariant functions A 
and B must be specified. To do this we assume that the conversion process proceeds 
via the mechanisms of Fig. 3, and consider as an example the transition amplitude for 
the ujN — » rjN process. Using the usual Feynman rules and vertices of the Lagrangian, 
Eqn. |l|, one writes, 

T^tjp = -MPi){9^NN*g V NN*l >1 l 5 [G R (s 13 ) + G R (ui 3 )) + 
guNN^Wl^GNisn) - 7 5 (^ - ta^k u )G N (u 13 )}}u(p 2 ) , (17) 

where Si 3 = (k + Pa) 2 , U\ 3 = (p 2 — k) 2 and Gjy denotes a nucleon propagator 

G N = 1- . (18) 

For an isobar the mass in Eqn. [18| is replaced by a mass operator, i.e., 

= 2 — n — , -p /o • ( 19 ) 
p- M R + iTr/2 

Generally speaking, Mr and Tr are the real and imaginary parts of the isobar self 
energy, quantities representing the mass and width of the resonance. Both of these 
are functions of the Mandelstam variable s = p 2 . It is straightforward to show that 
at resonance, where yfs ~ Mr, Eqn. |Tj] simplifies to the usual Breit-Wigner form, 

Gr ^ = M 7T T /9 • (20) 

Mr - Vs - iYr/2 
For future use we define two related quantities, 

A N (x) = i/(x — M 2 + ie) , (21) 

and 

A R {x) = —ij (Mr -V^~ ^ fi /2) . (22) 

Altogether there are four terms in Eqn. [T?| corresponding to s and u-channels with 
N* (diagrams 3a-3b) and with nucleon (diagrams 3c-3d) poles. By comparing Eqn. 
|17| with Eqn. [TT| one obtains, 

Aj.I = 9^NN* QriNN* G R + ig^NN— Xy , (23) 



Au,2 = -ZiMg^NN^^-GN 

m r , 



'if 



A w5 = i2Mg u , NN ^ ] ^-A N (si 3 ) , 



Similar calculations for other meson exchanges give the following: 

A n ■ frjNN v 
A-a — QaNN* QnNN* ^ R ~ ^QaNN -^S , 

B a = i 2Mg aNN ^G N , 

A rt -> frjNN v 
A-5 — g&NN*9r l NN* ( ^R'T ~ 1Q5NN A-S T 1 

B s = i2Mg SNN ^-^-G N T , 

A-,, = g^NN* QrjNN* GrT , 
D frjNN fnNN v -> 

±s n = Apr , 

Arj — 9rjNN*GR , 
f 2 

D JrjNN v 
&rj — 5— Ap , 

m r\ 

a rt -> , • fr\NN v -> 
A p,l — g P NN*grjNN*trRT + 'igpNN AyT , 



A Pi2 = -2iMg pNN ^Y v r , 

„ ■ frjNN K _ 

A P , 3 - -*9 P NN — ^Z V T , 
A p ^ = —2ig pNN — — kG^t , 



A pfi = %2Mg pNN 



^(l + 2 K )A N (s 13 y, 



A P:6 = -2ig pNN ^ vNN kA n (s 13 )t . 



In Eqns. - 40 we use the notations, 



G R = A R (s 13 ) + A R (u 13 ) , 
G N = A N (s 13 ) + A N (u 13 ) , 

X s = (2p 2 k + m 2 v - 2M 2 )A N (s 13 ) + (mj - 2 Pl k - 2M 2 )A N (% 
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X P = (2p 2 k + m 2 + 4M 2 )A^( Sl3 ) + (m 2 - 2 Pl k + 4M 2 )A^( Ml3 ) , (44) 

X v = (2p 2 k + m 2 )A N {s 13 ) + (2 Pl k + m 2 )A N {u 13 ) , (45) 
Y v = (1 + k + ^(^k + m 2 v ))A N (s 13 ) + 

(1 + « + 4^( 2 Pi fc - ^))A Ar (w 13 ) , (46) 

Z v = (2p 2 ■ k + m 2 TI )A N (s 13 ) - (2px • fc - m 2 )A w ( Ml3 ) . (47) 



2.4 The model parameters 

In the calculations to be presented below all of the meson-nucleon-nucleon couplings 
are taken from Machleidt et al.[[^|; their relativistic OBEP set with pseudovector 
coupling for the pseudoscalar mesons. For the sake of completeness these are included 
in Table |l[ As in Ref . J29| the source form factors are taken to be, 



Fa — 9oNNfa , F^ — ic/nN N fir , F v — ig v NNJr) ) F$ — ggNNfs 
Fp^ = g P NNfp , F^ = ^j^g p NNfp , Ffp = g u NNfu> , 



with the parameterization, 

M? 2 ) = • (48) 

The 7r and rj couplings to the N* (1535 MeV) are deduced from the partial widths 
of the N* to decay into 7rN and r/N; their values are rather well accepted (see Refs. 
[12-16]). The g P NN* is deduced from 7iV — > N* data and vector dominance [fP^] . For 



the other mesons there exits no direct information which can be used to determine 



their coupling to the N*. In Ref. [14] the u coupling is evaluated from the relation, 



g^NN* g-KNN* /^gx 



9uNN g-KNN 

Assuming SU(3) flavor symmetry it can be shown also that, 

g^NN* g P NN* 



(50) 



9wNN gpNNyl + K p ) 

Taking the appropriate constants from Table p] one obtains g^NN* ~ 1 from both of 
these expressions. This is a factor of two higher than that suggested by Vetter et 
al.[Qj. Similarly, from using effective triangle diagrams (see Appendix A) the g a NN* 
and gsNN* are also related to the pion constants through an expression similar to 
Eqn. gjj 
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3 Amplitudes and cross sections 

In this section we write S-wave production amplitudes for the pp — > pprj and pn — > pnr] 
reactions in a simple form which is suitable for numerical calculations. To this aim 
let us call, 

n, = p? , (51) 
Je~Tm 

where pj and Ej are three-momentum and total energy of the j-th nucleon. For the 
incoming particles in the CM system, IIx = —U 2 = II. The total energy square is 
s — (Pi + V2) 2 ■ Then the energy available in the CM system is, 

Q = y/lj - 2M - m v . (52) 

We shall calculate the amplitudes and cross sections as functions of Q. 

There are two amplitudes one isovector Mn and one isoscalar M 00 which deter- 
mine completely the cross sections for the NN — > Xn reactions at rest. These two 
amplitudes correspond to the 33 P — > 31 S and 11 Pi — > 31 S 1 transitions in the two- 
nucleon system. Only M n contributes to the rate of the pp — > pprj reaction. Likewise 
only M o contributes to the np — > drj reaction, but both amplitudes contribute to the 
rate of the np — > npn reaction. Following the discussion in subsection 2.2 we write 
these amplitudes as, 

M u = M n + M v + M a + M s + M p + , (53) 
Moo = -3M W + M v + M a - 3M 5 - 3M P + M w , (54) 

where Mb stands for a partial exchange amplitude of a meson B. The factor of —3 
is due to isospin. By evaluating matrix elements of the expressions, Eqns. 14-15, 
between the allowed initial and final spin states and substituting for the appropriate 
invariant functions listed in subsection 2.3, we obtain, 

M a = iG a NN(gaNN*9riNN*R + daNNdriNN^s) ? (55) 
Ms = iGsNN (95NN*9riNN*R + QSNN 9'qNN ^ 5 ) , (56) 
M n = iG n NN(gnNN*griNN*R + 9irNN 9riNN ^ p ) , (57) 

M v = iG vNN (gl lNN ^R + g^ NN T, P ] , (58) 
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+ 2v p )R + g pNN g vNN (£W + 2^)} , (59) 
± 2v u )R + guNNgnNNpW + 2E<?>)] . (60) 



In Eqns. 55-60 we have used the notation, 
„ E + M 1 

GbNN — gBNN 

R 



M M(m v + Q) + ml 
1 

+ 



f a (-M[m v + Q})U 
1 



M R - M -m v -Q -iT/2 M R - M + m v + Q - iT/2 



S = — I 1 - 5m ^ 
5 M I 2M- m r 



2M 



2(E + M) 

i + — ^n-n 



1 



2M + m v 



(1) _ 

' ~ 4M 2 



(i + «n • n) 



2 - a ^ m " 



E + M 



M 



n-n i 



2 ' 2M + m r 



Km; 



4M(2M + m, 



2M + m v J 4M(2M + m v ) 



m r 



2 - 



"AM 2 \~ 2M + m Vi 

sL 2) = o, 

w y = 2 + ( — ) n • n 



Vy = 1 + Ky 



2M 

E + M 
2M 



(61) 
(62) 
(63) 
(64) 

(65) 

(66) 

(67) 

(68) 
(69) 

(70) 
(71) 



The — and + signs in the expressions for M p and M w refer to the 33 Pq — *> 31 So and 
11 Pi — > 13 Si transitions, respectively. The expressions for other exchange amplitudes 
are the same for both of these transitions. Because of kinematic, the resonance terms 
in Eqns. 55-60 exceed by far any of the background terms. Following the definitions 
of w p and v p , the p exchange amplitude, Eqn. 59, has opposite signs for the two 
transitions but, because of the negative isospin factor the p contributions in Mn 
and M 00 have equal signs. Also note that for k u = 0, w u} =2 and v u =l causing the 
first term in Eqn. |60] to vanish in the case of a 33 Pq — > 3l So transition. Thus, only 
small background terms may contribute to the u exchange amplitude, and practically, 
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the calculated cross section for the pp — > ppt] reaction is expected to be free from 
ambiguities due to the u;NN* coupling. 

Finally a general formula for the invariant cross section for the NN — > NNr] 
reaction is written as 

°tt = }1 n ll 2 I d£l x dQ 2 ^X 1 %,mls2s) X 1/2 (s 23 ,ml,ml) \ M TT Z \ 2 , 
128 p*s 6 ' z J 

(72) 

where Mtt represents the amplitudes Mn and M 00 of Eqns. 53-54; the indices 1, 
2 and 3 label the outgoing particles, y/s and 1/523 are the total energy and partial 
energy of particles 2 and 3, p* is the center of mass (CM) momentum of the incoming 
proton, and Z is the FSI correction factor to be specified below. The A is the usual 



triangle function defined as [|33 



X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz . (73) 

4 Interference terms 

As indicated in the introduction, major uncertainties in previous model calculations 
arise from the unknown relative phases of different meson contributions. Particularly, 
interference terms between prominent p and 7r contributions influence strongly the 
scale of the cross section. The aim of the present section is to show that in the present 
model and in the limit where r^-meson production proceeds through N* (1535 MeV) 
excitations only, interference terms between the overall contribution from scalar and 
pseudoscalar meson exchanges and that of vector meson exchanges cancel out. To see 
this let us consider the contributions to the production amplitude from N* pole terms 
only. By substituting the various exchange contributions, Eqns. 55-60, into Eqn. 2 
one obtains, 

M™ N ^ PNN (pi,p 2 ;P3,P4, k) = g V NN*R 

[i(G n NNgTTNN* + G v NNgr)NN* + G a NNgaNN* + GfiNNgSNN* ) + 

GpNNQpNN* {w p ± 2v p ) + G pNN g pNN * (w u ± 2v w )} . (74) 
We notice that all of the quantities G bn n > g bn n* j wy and vy are real so that 

M (m) M (in)t = \ grjNm R\2 
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{[GttNN9wNN* + G-qNNg-qNN* + GaNNQaNN* + GsNNQSNN* ] ^ + 
[GpNNdpNN* {w p ± 2f p ) + GuNNQujNN* (™u ± 2« u )]" } . (75) 

We may thus conclude that in a model where isobar mechanism dominates, in prac- 
tice, the 7r and p exchange amplitudes add incoherently and their unknown relative 
phases should introduce no ambiguities in the model predictions. There remain of 
course ambiguities due to interference amongst the different vector meson exchanges 
or amongst the various scalar and pseudoscalar meson exchanges to be studied below 
in more details. 



5 Predictions 

We now apply the model presented in the previous sections to calculate the energy 
integrated cross section for the pp — > ppr] and pn — > pnr) reactions. We consider first 
the relative importance of various exchange contributions. To this aim we draw in 
Figs. 4-5 the partial exchange amplitudes, Eqns. 55-60. Most important are the 
p and 77 exchanges. Other contributions are small but might be influential through 
interference. Taking the relative phases of the different exchange amplitudes to be 
ones of signs, there are altogether 32 solutions, one for each of the different sign 
combinations. All of these solutions yield cross sections having identical energy de- 
pendence but vary in scale. The primary production amplitudes for the pp — > ppr) 
and pn — > pun reactions calculated with all of the relative phases chosen to be +1 are 
shown in Fig.6. As we shall see below this phase combination, hereinafter referred 
to as the standard solution, explains best the cross section data at energies close to 
threshold for the pp — > ppr], pn — > pnr). The cross sections corresponding to these 
amplitudes are drawn as solid lines in Figs. 7-8. The other curves in these figures 
represent predictions with other phase combinations, representing the lowest (dashed 
curve) and highest (dotted curve) cross sections obtained. At most the scale of the 
calculated cross sections varies by a factor of 2 — 3, as compared to a factor of about 20 
reported for other models[|l(J]. To examine these ambiguities in the model predictions 



let us consider the standard solution in some details. Consider first the pp — > ppr) 
reaction. Here the individual contributions from all of the scalar and pseudoscalar 
meson exchanges have equal signs and therefore should add constructively. In Fig. 
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9 we show partial cross sections accounting for scalar and pseudoscalar exchanges 
(dashed curve) and vector exchanges (dash-dotted curve) separately. The solid curve 
(labeled a + ) gives the total cross section for the standard solution. Although the size 
of these partial cross sections are comparable, reversing the relative phases for both 
of the p and uj to have the opposite signs hardly influences the predictions for the to- 
tal cross section (curve labeled u~). Indeed, interference between vector mesons and 
scalar or pseudoscalar exchanges involves background terms only and consequently 
are very small. Also, as we have indicated above the u exchange is very weak in 
this case because N* pole terms sum to zero so that the main source of ambiguities 
in the model predictions is due to interference amongst the scalar and pseudoscalar 
contributions. 

We obtain similar results for the pn — > pnrj reactions. In this case contributions 
from both of the 33 Pq — ► 31 Sq and 11 Pi — > 13 Si transitions are possible and the cross 
section becomes an incoherent sum, 

[cn + o"oo] • (76) 

The calculated cross sections are shown in Fig. 10. Again the phase between the over- 
all contribution from vector meson exchanges with respect to the overall contribution 
from scalar and pseudoscalar does not affect much the calculated total cross section. 
In this case the process is dominated by a p meson exchange. The contributions from 
uj exchanges through N* excitations do not cancel out (see Eqn. 60), giving rise to 
a relatively large uj exchange contribution, so that interference with a dominant p 
exchange amplitude determines, to a large extent, the scale of the calculated cross 
section. As for the pp — > pprj reaction, in the standard solution which explains data 
the p and uj add destructively. The isovector 7r and 5 contributions add destructively 
to those of the isoscalar rj and a. 

6 FSI and comparison with data 

At energies close to threshold FSI corrections influence the energy dependence and 
scale of the calculated cross section. To allow comparison with data to be made, we 



introduce FSI in an approximate way following the procedure of Refs. [16, [34 . 
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The S-wave transition amplitude for a three body reaction, say NN — > NNr], is 
assumed to factorize in the form 

T if « Ml]Tf f , (77) 

where Mfj stands for the primary production amplitude and T| j- represents the FSI 
correction factor. We identify this factor with the on mass-shell elastic scattering 
amplitude for the NNr] — > NNr] (three particles in to three particles out transition). 
The validity of this approximation was discussed at length previously jPJ]. Here we 
recall only that, Tft is a solution of the Faddeev's equations and has the structure 
of the Faddeev's decomposition of the t-matrix for 3 — > 3 transition. It is estimated 
using on mass shell two-body elastic scattering amplitudes, t v ^ and tjVTv, in a three- 
body space where the third particle being a spectator. We calculate these from 
S-wave phase shifts. The S-wave pp phase shift is calculated using the modified Cini- 
Fubini-Stanghellini ( CFS ) formula with a proton-proton scattering length a pp = 
— 7.82/m and an effective range r pp = 2.7 fm. Similarly, the np S-wave singlet and 
triplet phase shifts are calculated using the so called mixed effective range expansion 
with CFS shape corrections. The scattering lengths and effective ranges are taken to 
be at = 5.413 fm and r t = 1.735 fm for triplet scattering and a s = —23.715 fm and 
r s = 2.73 fm for the singlet. The values of the CFS shape correction parameters are 
taken to be p\ = 0.1147 fm 3 and p<i = 3.861 fm 2 . The t v w is taken from Ref. to 
be 

[t(rjN -+ rjN)]- 1 = 1/a + rk 2 /2 - ik v , (78) 
with ?7p scattering length and effective range 

a = (0.476 + 0.276i)/m, r = (-3.16 - 0.13i)/m . (79) 

Albeit, FSI corrections modify the energy dependence of the cross section at low 
energy region (Q below 15 MeV) and should not influence much predictions at higher 
energies. 

In Figs. 11-12 we draw our predictions for the energy integrated cross sections 
for the pp — > pprj and pn — > put] reactions with FSI corrections included, along with 
data taken from Refs. [2-6]. The solid curve gives the standard solution, i.e., with 
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the relative phases of the various exchange amplitudes taken to be +1. The model 
explains fairly well the cross section data for both of these processes. Calculations 
which include resonance terms only and without background terms are drawn as 
dashed curves. Clearly, the n production is dominated by the N*(1535 MeV) isobar 
excitations (diagram 2a). The effects of the background terms (diagrams 2b-2d) 
are bound to 10% or less in the near threshold energy region considered here. FSI 
corrections for the np — > npn reaction are performed separately for an and <Too using 
the appropriate singlet and triplet np S-wave phase shifts. The model reproduces the 
cross section data for all of the pp — > ppn, pn — > pnn and pn — *> dn reactions. The 



results for the latter process will be reported elsewhere p6[ . 

As indicated in the introduction, there exits no direct evidence for the coupling 
of the a, 5 and u mesons to the N* (1535 MeV) isobar. It is interesting to see 
whether contributions from these mesons are required to explain data. To this aim, 
calculations were repeated with : (i) g^NN* = 0, (ii) g^NN* = g^NN* = 0, (iii) 
g^NN* = gaNN* = gsNN* = 0. The predictions with g^NN* = are shown in Figs. 
11-12 as dot-dashed curves. They explain equally well the data for the pp — > ppn 
reaction, and for the pn — > pnn reaction the agreement is even slightly better. The 
effects of the a and 5 are negligibly small and their contributions seem not to be 
needed either (see Figs. 12-13). 

7 Summary and Discussion 

In this paper we have presented a covariant OBE model for n production in NN 
collisions. The model includes nucleon and N* (1535 MeV) isobar excitations in 
the intermediate states and contributions from all of the ir, n, a, p, uo and 5 meson 
exchanges. Our starting point was a Lagrangian which describes NN scattering, in a 
consistent manner, in terms of these same meson exchanges. 

As in other similar calculations, the model predictions are subject to uncertain- 
ties due to the unknown relative phases of the various exchange amplitudes, but 
these are now limited to only a factor of 2-3 in the calculated cross sections. The 
present formalism allows a thorough discussion of interference between the various 
contributions. We have shown that interference effects between vector and scalar or 
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pseudoscalar meson exchanges are negligibly small. Bearing in mind that the most 
prominent contributions are those due to the 7r and p exchanges this feature of the 
model renders predictions to be considerably more reliable than our previous model 
calculations |L6|]. Assuming that the relative phases of all of the various exchange 
amplitudes to be +1, the model reproduces the cross section data, both the scale and 
energy dependence, for the pp — > pprj and pn — > pnrj reactions. As indicated above 



the model also reproduces the cross section data for the np — > drj reaction [[36 



The u, a and 5 contributions are small and with the present accuracies of the 
measurements seem unnecessary to explain the data. This is not unwarranted prop- 
erty of the model since there is no direct evidence for the coupling of these mesons to 
the N* (1535 MeV) isobar. Finally, with the p coupling lowered to g P NN* = 0.615, as 
low as suggested by Vetter et al.fll4]|, no solution is found which explains data simul- 
taneously for all of the AW — > Xrj reactions. It would be interesting to determine 
this coupling anew. 
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8 Appendix A 

In this Appendix we describe the derivation of the expressions used to calculate g a NN* 
and g$NN* ( see caption of Table 1). Consider the effective triangle diagram of Fig. 
14 and take the crmi vertex to have the general form, 

V a7T7T = g a7T7T v a7T7T (k, q-k,q) . (80) 

For a Lagrangian density, L int = g ulT1T MaTi ■ n, the function v an7T (k, q — k,q) becomes, 

v (77T7T (k,q-k,q)=Mf l {k 2 ,q 2 ) , (81) 
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where fi(k , q ) is a scalar function. Taking the Lagrangian density to be, L 
(g a7T7T /M)(adif) ■ (07?), yields, 

Vann(k, q - k, q) = —(q - k)qf 2 (k 2 , q 2 ) . (82) 

As will be demonstrated below the actual form of f CT7r7r does not influence the final 
results for the calculated coupling constants. Suppose now that one of the nucleon is 
on mass shelly. e., p 2 = M 2 , then using the loop diagram Fig. 14 the aNN coupling 
is 



37] 



gaNNVaNN(k,p - k,p) = -% ^ T 4 J d A q^v a7T7T (k, q-k,q) 

q 2 - ml)[(q - k) 2 - ml][(p - q) 2 - M 2 }}- 1 . (83) 

By using the standard Feynman parameterization and shift of the integration variable 
one obtains, 

g a NNV a NN(k,p - k,p) = g an7T gl NN [p'X(k,p - k,p) + ftY^kiP - k,p)] , (84) 
where, 

X{k,p-k,p) = 2fdx 1 dx a dx 3 x s 8{l-x 1 -x 2 -x i ))I q , (85) 
Y(k,p — k,p) = 2 J dxidx 2 dxsX25(l — x\ — x 2 — x 3 ))I q , (86) 

and I q is the usual integral over the Euclidean four- vector g[31] , 

I q = t^u J d^qv^k, q + px 3 + k[x 2 -l],q + px 3 + kx 2 ) 

(q 2 — (q + px 3 + kx 2 ) 2 + k 2 x 2 + m^{x\ + x 3 ) . (87) 
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To evaluate the g a NN* coupling constant we have to replace one of the two irNN 
vertices with a nNN* vertex. This leads to, 



gaNN*V aNN * (k, p-k,p) = g^girNNQirNN* \j>X{k, p - k,p) + $Y(k, p ~ k, p)} . (88) 

After some algebra, 

9aNNVaNN{k,p - k,p) = g U7T7T gl NN J , (89) 

g<rNN*V aNN * (k, p- k,p) = gawTrg-KNNg-KNN* J , (90) 



with, 



J 



M 2 X 2 + k 2 Y 2 + 2pkXY 1/2 . (91) 



From Eqns. 89 and 90, one sees that the functions v a NN and v^nn* are identical up 
to a scale factor, which we may absorb into the constant g a NN*- Thus, Eqns. 89 and 
90 are equivalent to, 

gaNN* _ g-KNN* ^g 2 ) 
gaNN g-KNN 

Note that taking the a coupling to be ig a NN*[N*^ b N + iV7 5 iV*], as in Eqn. [I], lead 
to a real g a NN*- Otherwise with this coupling taken as g a NN* [N*rfN — Nj 5 N*}, one 
obtains a pure imaginary g a NN*- 

Similar calculations can be made for the gsNN* leading to a similar expression, 

viz, 

' f ' " ' ' ' (93) 



gSNN* _ l 
gSNN 2 

9 Appendix B 



g-KNN* g-qNN* 



gKNN g-qNN _ 



In what follows we write explicit expressions for the vectors and tensors of Eqns. 
14-16. Let i,j label the i-th and j-th nucleon (see Fig. 1). We define, 

^liLVri,, (94) 

Eij = U { + Uj , (95) 

Ay = Ili - n, , (96) 

Xij = Ik x LL, . (97) 
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Then for a scalar meson exchange, 

S = tF s (q 2 )G s (q 2 )T^ , (98) 

S[ 3 = Fs(q 2 )G s (q 2 )TPx[3 , (99) 

S l 24 = F s (q 2 )Gs(q 2 )<j>i 3 TP , (100) 

S kl = tF s (q 2 )G s (q 2 )x k 13 Tp i , (101) 

with, 

TP = ±(k + tj)-Xi3B s , (102) 

Tf = a s a 1 24 + ±B s [(k + q ) • n 4 n 2 + (k + q ) • n 2 ni + 

2 - 4 (k + q)'-(A;o + ?o)Sy • (103) 
For a pseudoscalar exchange, 

P< 3 = G P (q 2 )F P (q 2 )T^A[ 3 , 

P fc/ = zF P (g 2 )Gp(g 2 )At 3 Tf 1 , (104) 



with, 



T p 1} = Apfa + ^ 5 [0+(A;o + g ) + (k + q) • £ 24 ] , 

T p 2)l = A PX l 24 - l -B s [(k + g ) X24 + (k + q) x A 24 ] . (105) 



For a vector meson exchange, 



where, 



V = tG v (q 2 )C? 3 X 2 ^ , 

^i3 = -Gv{q 2 )Dilx 2A)lJi , 

V 2 \ = G v (q 2 )C^Y 24 fM , 

f to = iG v (q 2 )D^Y 2 ^ , (106) 

Ci° 3 = 4 1 W4 + 4 2) q-Ai3 , (107) 
cs = f^s™ + 4 2) ( g °Ar 3 - n^q • n 3 - n™ q • no , (ios) 
= -xL4 1} + q x ^ 3 4 2) , (109) 
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D rnk = e m« [(A r 3F a) + g 0^ + 0+ g^f)^ - q r (U^ + U{U^)F^] , (110) 

X 2 ° 4 = k • X24 a° 6 , (111) 
X% = k • X 24 < 6 - e mrs n^n^x - e mrs FA^ 2 , (112) 

1m'' = («°5 + k°A 2 )A l 24 - a%[k°X l 24 - <f) 24 k l + U l 2 U 4 • k - U l 4 U 2 • k] 

- aj: 1 24 + a 2 [A;V 2 - 4 + n 2 n 4 • k + n 4 n 2 • k] , 

Ym 1 = « + k m A 2 )A l 24 - < 6 (£; E 24 - 24 A; Z - Tl l 2 Tl 4 • k - Tl l 4 Tl 2 • k) 

- A 1 (<j) 24 5 ml + n™n 2 + n™n 4 ) - A 2 [k • A 24 5 ml - A™k l 

- k°((j) 24 5 ml + n™n 2 + n™n 4 )] . (113) 
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Table 1: The model parameters. The constants gsNNi cut off parameters (A#) and 
meson masses (mg) are taken from Machleidt et al. ||29[| ; their relativistic OBEP solu- 
tion with pseudoscalar coupling for tt and r\ and with f P NN Iq p nn = 6.1, f w NN / QuNN = 
. 





m B (MeV) 


9bnn/^ 


A B (MeV) 


Qbnn* 


71 


138 


14.6 


1300 


0.8 


a 


550 


8.03 


1800 


0.5 b ) 


7] 


547.3°) 


3.0 


1500 


2.2 


P 


769 


0.95 


1300 


1.66 


uo 


783 


20.0 


1500 


0.94 6 ) 


5 


983 


5.07 


1500 


1.48 c ) 



a ) From Ref.|| . 

fc ) From g a NN*/g<TNN = 9nNN* / QivNN 

c ) From gsNN* /gsNN = \ [g^NN* / g-wNN + g-qNN* i 'svvtv] 




Figure 1: The primary production mechanism for the NN — > NNt] reaction. A bo- 
son B created on nucleon 1 (momentum pi), is converted into an r\ meson (momentum 
k) on nucleon 2 (momentum p-i). 
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Figure 2: Feynman diagrams for r] meson production in NN collisions, (a) s-channel 
N* pole contribution; a boson B (dashed line) formed on nucleon 1 (momentum pi) 
is absorbed on nucleon 2 (momentum P2) which is excited to an isobar state with 
momentum p = k + p 4 which then decays into a nucleon 4 (momentum p 4 ) and an r\ 
meson (momentum k). (b) u-channel N* pole; pre-emission counter part of diagram 
(a) where an 77 meson is emitted before collision, (c) s-channel N pole; as in (a) 
with a nucleon in the intermediate state, (d) u-channel N pole contribution; pre- 
emission counter part of (c). There are four more graphs in which the exchanged 
boson is formed on nucleon 2. These are not displayed here but their contributions 
are included in the calcultions. 
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Figure 3: Feynman diagrams for the conversion process BN — > r]N (a) s-channel N* 
pole contribution, (b) u-channel N* pole; pre-emission counter part of diagram (a), 
(c) s-channel N pole; as in (a) with a nucleon in the intermediate state, (d) u-channel 
N pole contribution; pre-emission counter part of (c). 
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Figure 4: Partial exchange amplitude for the pp — > pprj reaction. The uj contributions 
from s and u-channels with N* poles cancel out and is very weak (see text). 
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Figure 5: Partial exchange amplitude for the pn — *> pnt] reaction. 
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Figure 6: Primary production amplitude for the pp — > ppr\ and pn — > pnry reactions 
with all of the relative phases taken to be +1 (our standard solution). 
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Figure 7: Predictions for the pp — > pp?7 cross sections near threshold with differ- 
ent phase combinations. The solid line displays the standard solution. The lowest 
(dashed curve) and the highest (dash-dotted curve) cross sections are obtained with 
the n, T], a, 5, p and u> relative phases being — h + + ++ and 1 — , respectively. 
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Figure 8: Predictions for the pn — > pni] cross sections near threshold with differ- 
ent phase combinations. The solid line displays the standard solution. The lowest 
(dashed curve) and the highest (dash-dotted curve) cross sections are obtained with 
the 7i, r], a, 5, p and u> relative phases being + + H h+ and — h H h, respectively. 
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Figure 9: Interference between vector and scalar-pseudoscalar meson contributions. 
The curves labeled by oy, cp+s, represent partial cross sections with vector meson 
exchanges and with scalar plus pseudoscalar exchanges only. The other lines are total 
cross sections with all the relative phases taken as +1 and with the same relative 
phases for the scalars and pseudoscalars except for M p and M u having the oposite 
sign (a-). 
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Figure 10: Interference between vector and scalar-pseudoscalar meson contributions 
for the pn — > pnr) reaction. See captions of Fig. 9. 
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Figure 11: Energy integrated cross sections for the pp — > ppr] reaction versus 
the energy available in the CM system. All curves include FSI corrections via the 
approximation Eqn. 77. Predictions are given for the standard solution with both 
resonance and background terms included (solid line) and with the resonance term 
only (dashed curve). The dot-dashed curve (which in this case override the solid line) 
gives the standard solutions as obtained with g w = 0. The data shown are taken from 
Refs. 2-4,7. 
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Figure 12: Energy integrated cross section for the pn — > pnr/ reaction. Data point 
are taken from Ref. 7. See caption of Fig. 11. 
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Figure 13: Effects of u, a, 5 exchanges. Predictions for the pp — > ppr] cross section are 
given for the standard solution with (i) q^nn* = (solid line), (ii) g^NN* = QoNN* = 
(dashed curve) and (iii) Qunn* — 9<tNN* — Qsnn* =0 (dash-dotted curve). 
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Figure 14: Effects of u, a, 5 exchanges. Predictions for the pn — > pnrj cross section are 
given for the standard solution with (i) q^nn* = (solid line), (ii) g^NN* = QoNN* = 
(dashed curve) and (iii) Qunn* — 9<tNN* — Qsnn* =0 (dash-dotted curve). 



39 



a 



I k 



V 



\ q - k 
\ 



JT 



P - k 



N 



p - q 



N (N*) 



Figure 15: An effective triangle diagram used to derive an expression for g a NN* 
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